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ABSTRACT 

We present a feedback mechanism for supermassive black holes and their host bulges that operates 
during epochs of radio-loud quasar activity. In the radio cores of relativistic quasar jets, internal 
shocks convert a fraction of ordered bulk kinetic energy into randomized relativistic ions, or in other 
words cosmic rays. By employing a phenomenologically-motivated jet model, we show that enough 
1-100 GeV cosmic rays escape the radio core into the host galaxy to break the Eddington limit in 
cosmic rays. As a result, hydrostatic balance is lost and a cosmic ray momentum-driven wind develops, 
expelling gas from the host galaxy and thus self-limiting the black hole and bulge growth. Although 
the interstellar cosmic ray power is much smaller than the quasar photon luminosity, cosmic rays 
provide a stronger feedback than UV photons, since they exchange momentum with the galactic gas 
much more efficiently. The amount of energy released into the host galaxy as cosmic rays, per unit 
of black hole rest mass energy, is independent of black hole mass. It follows that radio-loud jets 
should be more prevalent in relatively massive systems since they sit in galaxies with relatively deep 
gravitational potentials. Therefore, jet-powered cosmic ray feedback not only self-regulates the black 
hole and bulge growth, but also provides an explanation for the lack of radio-loud activity in relatively 
small galaxies. By employing basic known facts regarding the physical conditions in radio cores, we 
approximately reproduce both the slope and the normalization of the M, — M^, relation. 
Subject headings: cosmic rays — galaxies: formation — galaxies: fundamental parameters — galaxies: 
jets 



1. INTRODUCTION 

There are two classes of persistent sources at cos- 
mic distances: galaxies and quasars/ Active Galactic Nu- 
clei (AGN). Stars power galaxies while accretion onto 
and/or spin down of supermassive black holes power 
quasars. Until recently, galactic and quasar phenomena 
were thought to be separate on both observational and 
theoretical grounds. However, the discovery of the black 
hole mass - bulge stellar mass relation (M, — M^) in 
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and the black hole mass - stellar velocity dispersion rela- 
tion (M. -cr^; | Gebhardt et al. | 2000||Ferrarese fc Merritt 
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indicates that galactic and black hole activity are closely 
connected to one another. The natural implication is 
that the energy release resulting from the build-up of 
the black hole mass limits any further growth of both 
the stellar bulge and the black hole. The fact that the 
M, — relation holds for nearly four decades in black 
hole mass seems to suggest that a universal, self-similar 
or scale-free process is at work, which acts to self-regulate 
the ratio between black hole and bulge mass, irrespec- 
tive of their combined mass. Apparently, the only ques- 
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tion that remains is with regard to the exact physical 
mechanism responsible for black hole feedback and self- 
regulation. 

The Soltaiil(|1982| argument, along with the work of |Yu| 
|fc Trem ainc (2002), indicates that the mass of supermas- 
sive black holes is mostly accrued during an optically- 
luminous radiatively-efhcient "quasar phase." The en- 
ergy released during the accretion process, which is car- 
ried away primarily by photons and/or a "quasar wind," 
may couple to the interstellar medium of the host galaxy 
and eject it from the galactic gravitational potential (e.g., 
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In doing so, fuel for any further 
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galactic and quasar activity is removed and the mass of 
the black hole, as well as of the stellar bulge, is self- 
limited. 

The energy released during the accretion process may 
be carried away not solely in the form of photons. In 
the so called "radio-loud" (as opposed to "radio-quiet") 
objects, relativistic coUimated outflows, or "jets," put 
out a significant amount of energy in mechanical form. 
Although radio-loud phenomena are also ob served in ob- 



jects that are n ot actively accreting (e.g., Ho & Peng 



2001| Ho 20021, the radio jet is more likely to attect 
the evolution of the system when a significant amount 
of mass is being built up. As this work focuses on the 
self-regulation of black hole growth, which occurs at high 
accretion rates, our attention rests on objects that are 
both significantly accreting and radio-loud. 

The kinetic power of radio jets is dissipated in sub-pc 
scale "radio cores" and kpc to Mpc scale "radio lobes," 
with comparable amounts of energy dissipated at each 
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site. The radio-loud quasar phase could be responsible 
for black hole self-regulation if the energy release from 
the radio core, unlike that from the distant radio lobes, 
has the opportunity to couple to the interstellar medium 
of the host galaxy. 

In this work, we propose that black hole self-regulation 
may be mediated by 1-100 GeV protons, or in other 
words cosmic rays, produced in the jet core during phases 
of radio-loud quasar activity. We show that enough cos- 
mic rays escape the radio core into the host galaxy to 
power a cosmic ray-driven wind that ejects the interstel- 
lar gas, thus removing the fuel for further star formation 
and black hole accretion. 

Our arguments are organized as follows. In S|2]we lay 
out the basic idea behind our jet-powered cosmic ray 
feedback scenario by defining and recognizing the im- 
portance of the Eddington limit in cosmic rays. In ^|3] 
we briefly summarize the physical features of the stan- 
dard model of radio-loud AGN jets that are important 
for our feedback mechanism. With the jet model in hand, 
in ^|4] we determine the interstellar cosmic ray luminos- 
ity resulting from radio-loud quasar activity in terms of 
observable quantities. For readers uninterested in the 
details of our jet and radio core model, as well as its 
phenomenological underpinnings, we suggest skipping ^|3] 
and f|4] In J 5] we examine the consequences of our cosmic 
ray feedbacK scenario and make a critical comparison of 
the black hole self-regulation model presented here with 
other models of quasar feedback. We summarize our find- 
ings in 36] where we also present a brief review of well- 
known observations that lend support to our jet-powered 
self-regulation mechanism. In addition, we propose an 
observational test. 



efficiently couple to the gaseous component of the host 
galaxy. In other words, radio-loud kinetically-dominated 
epochs may represent another mode of AGN activity, 
other than optically-luminous radiatively-efhcient quasar 
phases, that could potentially lead to efficient black hole 
self-regulation. 

The photon luminosity of the sub-pc scale jet core re- 
sults from the cooling of relativistic electrons through a 
combination of synchrotron and inverse Compton emis- 
sion. By modeling the spe ctral energy d istribution of 
powerful radio-loud quasars, Celotti & GhiscUini (2008| 
infer that the emitting electrons must be accclcrEiEea, 
presumably by some magnetized Fermi mechanism, up 
to highly-relativistic energies, in excess of a few tens of 
GeV (in the jet comoving frame) . The details of the ac- 
celeration mechanism aside, it is reasonable to conclude 
as well that protons are randomized with energies up 
to ^ 1 — 10 GeV. If a significant fraction of these ran- 
domized relativistic ions, or in other words cosmic rays, 
escape the site of acceleration and diffuse into the inter- 
stellar medium of the host galaxy, they may provide the 
coupling between jet power and galactic gas required for 
black hole self-regulation during radio-loud phases. 

For example, in actively star-forming galaxies the gen- 
eration and subsequent diffusion of cosmic ray protons 
- a well-known by-product of core-collapse supernovae, 
and thus massive star-formation - may act as an agent 
of self-regulation, limiting the rate of star-formation and 
therefore the galaxy luminosity. The simplest way to un- 
derstand this is by defining the Eddington limit in cosmic 
rays 



2. THE EDDINGTON LIMIT IN COSMIC RAYS 

If supermassive black hole growth res ults primarily 
from radiatively-eflicient accretion, as the Soltan ( 1982 1 
argument suggests, then the ratio of the total energy re- 
leased during the accretion process to the binding energy 
of the galaxy's gaseous phase is given by 



AE, 

1^ 



M,/A'E 

10-3 



f-1 -2 
-1 Jg,-! "*,300 



,(3) 



galaxy system follows the 
Un. Here, e„^ — 0.1 e 



in case the black hole ■ 
M, — Mi, relation in eq. ^^^^ — 

is the radiative efficiency of the accretion ffow. The 
binding energy of the gas Eg ~ fgM^^a'^ is appropri- 
ate for a galaxy resembling an isothermal sphere whose 
gravitational force is mainly provided by stars and dark 
matter; fg = 0.1/g^_i is the galaxy's gas fraction and 
(T^ = 300 (T^_3oo kms~^ is the stellar velocity dispersion. 
Clearly, only a small amount of the photon energy re- 
leased during optically-bright quasar epochs is needed 
to couple to the galactic gas in order to eject it from 
the galaxy's gravitational potential, thus self-limiting the 
combined growth of the black hole and the stellar bulge. 

Again, during radio-loud phases a significant amount 
of the energy resulting from black hole accretion is re- 
leased in mechanical form as powerful relativistic jets. 
Even if the time-integrated kinetic output of the jet is 
AEj ^ AE,, eq. dsl) indicates that there may be ample 
energy to unbind the galactic gas, provided that the frac- 
tion of jet kinetic energy deposited in the radio core can 
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where L^^^ ^ is the Thomson Eddington limit for a galaxy 
with mass AI^,; we have pinned the ratio between the 
cosmic ray mean free path and the Thomson mean free 
path to the value A^^p^/A^, ~ 10~^, as in the case of 
the Milky Way. An interstellar cosmic ray luminosity 
of L^,jj ~ 10'**ergs~^ - equal to the corresponding cos- 
mic ray Eddington limit L^^^ ^-^ for the most massive 
galaxies - may result from the act of forming stars at 
a rate of ~ 10^ Mq yr~^ - corresponding to the bright- 
est star- forming galaxies. The formation of stars at a 
higher rate leads to the breaking of the cosmic ray Ed- 
dington limit and a cosmic ray-driven wind develops, 
removing the gaseous phase of the galaxy thus choking 
off star-formation. From this straightforward argument, 
it is quite possible that the luminosities of star-forming 
galaxies are capped at their observed values because they 
are Eddington - limite d in cosmic rays, as proposed by 



Socrates et al.| ( |2008 ). 

With respect to AGN, radio-loud objects display pow- 
erful jets with kinetic luminosity upwards of L, ~ 



lO'^^ergs ^ (e.g., Rawlings fc Saunders| |l991 1, compara- 



ble to the Thomson Eddington limit Z^^J^Tor the most 
massive black holes with M, ~ 10^ Mq. If the jet kinetic 
power is a fraction A^^^ ~ 1 of the black hole Thomson 
Eddington limit, we have 
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where M, = lO^M.^gM©. Let be the efhciency for 
converting the jet kinetic power Lj into interstellar cos- 
mic ray power L^j^ . A statement of momentum conser- 
vation and hydrostatic balance, the Eddington limit in 
cosmic rays defined in eq. Q indicates that a galactic 
cosmic ray-driven wind develops if 

^CR = <^CR -^.T ^ -^Edd.CR 1 (6) 

or if the fraction of jet kinetic luminosity ending up as 
interstellar cosmic rays attains a value of 
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for a black hole - galaxy system that follows the Af, — 
relation. 

The momentum requirement in eqs. ([6| or ([t]) is a 
necessary condition to initiate a galactic cosmic ray- 
driven outflow. However, in order to fully unbind the 
galaxy's gaseous phase, thus self-regulating the black 
hole and bulge growth, this epoch of super-Eddington 
cosmic ray activity should last long enough so that the 
time-integrated energy injected as cosmic rays into the 
interstellar medium is comparable to the binding en- 
ergy of the galactic gas. Interestingly, the value of the 
jet-cosmic ray conversion efficiency given in eq. ([7| is 
roughly equal to Eg//S.E, for the most massive galax- 
ies (see eq. pi), where Eg/AE, may be viewed as the 
minimum efliaency of a "generic" feedback mechanism 
for coupling the black hole energy release to the gaseous 
component of the host galaxy. In other words, the min- 
imum value of required to break the galaxy's mo- 
mentum balance, given in eq. ([t]), is approximately equal 
to the minimum value of the energy coupling efficiency 
Eg/AEt for the most massive systems.^ 

Clearly, an understanding and determination of the jet- 
cosmic ray efficiency parameter e^j^ is of fundamental im- 
portance with respect to quantifying whether or not jet- 
cosmic ray feedback is in fact the agent of self-regulation 
for a given black hole ~ galaxy system. In what follows, 
we take into account the physical properties of AGN ra- 
dio cores with the help of a well-established jet model. 
In doing so, we are able to determine the jet-cosmic ray 
efficiency parameter e^^^^ and discuss the effectiveness of 
our cosmic ray-driven feedback scenario. 

3. A STANDARD JET MODEL FOR RADIO-LOUD AGN 

Here we describe the physical properties of quasar jets, 
whose radio core is responsible for injecting cosmic rays 
into the interstellar medium of the host galaxy. The 
reasoning and details of the jet model outlined below 
are grounded by decades of observations of radio-loud 
quasars. We also discuss some of the observational con- 
straints that we use to determine the parameters of our 
model. 

3.1. Basics of Radio-Loud AGN Phenomena 

The spectral signature of radio-loud AGN results from 
dissipation of the kinetic energy of powerful relativis- 
tic jets in sub-pc scale "radio cores" and kpc to Mpc 

^ Strictly speaking, the minimum energy coupling efficiency for 
our jet-powered cosmic ray feedback would be ~ Eg/AEj, where 
Ai?j is the time-integrated kinetic output of the radio jet. We are 
implicitly assuming that AEj ~ AE, at the high-mass end. 



scale "radio lobes." The basic radio- loud phenomenology 
is lar gely understood in terms of a viewing-a ngle effect 
Antonucci 1993 Urry fc Padovani|1995 ). Emission 
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from the relativistic jet flow in the radio core is highly 
beamed due to relativistic aberration. The radiation seen 
along the jet axis is Doppler-boosted in both frequency 
and flux, when compared to the flow rest frame, and the 
opposite is true when the radio core is viewed edge-on. 

Core-dominated on-axis sources are collectively re- 
ferred to as "blazars" ~ a marriage between BL Lacs 
and flat-spectrum radio quasars (FSRQs). They dis- 
play extreme variability, broad (radio to gamma-ray) 
spectral energy distributions, super-luminal motion and 
super-Eddington fluxes. BL Lac spectra are typically 
featureless, whereas intense line emission is detectable 
in the optical-UV spectrum of FSRQs, presumably aris- 
ing from the accretion disk and its associated broad-line 
region (BLR). The continuum spectral energy distribu- 
tion of blazars shows two prominent bumps, with the 
high-energy bump that often dominates the bolometric 
power output. The low-energy bump, at infrared to UV 
frequencies, is conventionally attributed to synchrotron 
radiation, while the high-energy emission, extending up 
to gamma-ray energies, is thought to result from inverse 
Compton up-scattering of synchrotron photons (in BL 
Lacs) or external photons from the accretion disk and 
the BLR (in FSRQs). 

In objects where the jet is viewed edge-on, the so called 
"Fanaroff- Riley" galaxies, the radio emission mostly re- 
sults from the dissipation of jet power in radio lobes on 
inter-galactic or super-galactic scales. The emission from 
sub-pc scales is largely beamed away from the observer 
for these off-axis sources. In the framework of unification 



models (e.g., Antonucci 1993 Urry & Padovani 19951, 
FR I radio galaxies are identified as the parent popula- 
tion of BL Lacs, while FR II galaxies are the edge-on 
counterpart of FSRQs. 

In this work, all of our attention will be directed to- 
wards the most powerful objects, i.e., FSRQs and FR II 
galaxies, sources where radio- loud activity coexists with 
substantial black hole accretion. From here on, FSRQs 
and their parent population of FR II radio galaxies will 
be collectively referred to as "radio- loud quasars" . 



3.2. A Simple Jet Model 

We assume that the jet power or "luminosity" is 
primarily carried in mechanical form by matter expelled 
from the central engine at a rate Mj , such that 



TA-Ec^ 



(8) 



where F is the characteristic Lorentz factor of the flow. 
Under the assumption that the ultimate source of jet 
power is accretion onto the black hole at a rate M, 
(as opposed to extraction of the black hole spin), then 
Mj / M, ~ e^^j /r if the jet kinetic power is comparable to 
the accretion photon luminosity, ~ L^^^ ~ e^^^M^c^. 
It follows that for typical values, n amely e^,^^ ~ 0.1 and 
r - 10 (e.g., [Jorstad et al.|[2005| , we have MJM, ~ 
0.01. In other words, a significant fraction of the ac- 
cretion power must be transmitted to only ~ 1% of 
the mass. The mechanism of this seemingly improba- 
ble transfer of energy is assumed to be an unknown from 
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here on as we are primarily interested in its mechanical 
output. 

Powerful radio-loud quasars are variable, often display- 
ing radio blobs, or "knots," moving along the jet. Such 
behavior is commonly thought to result from the dissipa- 
tive interaction of shells of matter intermittently ejected 
from the central engine. If the Lorentz factors of the 
shells differ such that AF ~ F in the reference frame of 
the host galaxy (here, F may be thought of as the mean 
Lorentz factor of the jet flow), then the characteristic dis- 
tance R^-^^ from the black hole where dissipation takes 
place is 



comoving frame 



-R^ ~ 0.05 F^e-;;;! pc 



(9) 



where F = lOFi; R^ — GM,/c^ is the black hole gravi- 
tational radius and R^/e^^^c is the dynamical timescale 
of the accretion flow, which determines the typical time 
delay between subsequent shell ejections by the central 
engine. For the sake of simpUcity, we adopt a one-zone 
model in which the whole jet spectral energy distribution 
originates from a single region located at R^^^^ ■ 

We assume that the jet is b eamed into a cone with 
half-opening ang le 6* ~ 0.1 (e.g., |Jorstadet al.|2005[ ). For 
typical values F '--^ 10 and 9 ~ 0.1, in the jet conioving 
frame the dissipation region resembles a spherical blob, 
since its comoving longitudinal length ~ R^iBs/^ is com- 
parable to its transverse size ~ ^dian^- For a bi-conical 
mass-dominated jet. 
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where R is the distance from the black hole and p is 
the jet mass density at radius R measured in the host- 
galaxy frame. If the jet kinetic power is a fraction Aj^^^ of 
the black hole's photon Eddington limit, the density p^.^^ 
at the dissipation scale R^^^^ and the Thomson optical 
depth Tjj^^ down to R^^^^ niay be written in terms of basic 
physical parameters as 
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where p^.^ — p^.^^/F is the jet mass density at R^^^^ mea- 
sured in tlie comoving frame. If the two colliding shells 
have equal rest mass and differ by AF in Lorentz fac- 
tor, the fraction of bulk kinetic energy converted into 
heat is e., ~ 1/8 (AF/F)^ to leading order in AF/F (e.g., 



Kobayashi et al. 
tor shells with Ai' 



1997|. We take e^^^ ~ 0.1 as a benchmark 
F. It follows that the average co- 
moving Lorentz factor of shocked protons is 7p ~ 1 -|- e^^^ , 
corresponding to a mean energy ~ F 7p WpC^ 10 GeV 
in the frame of the host galaxy. If a sufficient number of 
these randomized relativistic ions, or cosmic rays, escape 
the jet into the galaxy's interstellar medium, they may 
provide the coupling between jet power and galactic gas 
required for our feedback model. 

Since radio cores are inarguably powered by syn- 
chrotron emission, our jet model cannot be com- 
plete without an estimate of the magnetic field. We 
parametrize the magnetic energy density in the comov- 
ing frame as a fraction '--^ 0.1 of the thermal energy 
density U'^^: 

U' . (15) 

B B th ^ ' 

It is reasonable to suspect that the magnetic field is 
tangled and inhomogenous within the dissipation region, 
more so since we assume the jet electromagnetic flux to 
be sub-dominant with respect to the kinetic flux. In the 
comoving frame, the mean magnetic field strength is 

B' ~ 2.7 , ^-3 e-Hl,l'^ G , (16) 

B,-l th,-l rad,-l Edd i -i J V / 



where e,, = 0.1 e . and e„ = 0.1 e , . 

th th,-l B B,-l 

Now, we possess all of the necessary physical ingredi- 
ents for our feedback model. However, before we assess 
whether or not a sufficient number of cosmic ray pro- 
tons are able to leak out of the radio core, potentially 
leading to the disruption of the entire host galaxy, we 
spend some time - for the sake of completeness ~ on the 
observational motivations that went into our jet model. 
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0.02-e^^^_^A,,,Fr^ri' , (12) 



where 9 — 0.1 0_i and k^s — 0.4cm^g"^ is the electron 
scattering opacity; Up and rig are the jet number den- 
sities of protons and leptons (including both positrons 
and electrons) in the host-galaxy frame. In eq. ( 12 1 we 
assume that the jet kinetic power is mostly carried by 
protons - in |3.3.2| we estimate that pairs outnumber 
ions by o nly a factor of ten. For rie/rip ~ 10, it follows 
from eq. ( 12 1 that the dissipation region is optically thin. 



In the dissipation region, internal shocks resulting from 
shell-shell collisions convert a fraction e^^ of the orderly 
mechanical energy of the jet flow into heat. It follows 
that the randomized thermal power generated in the dis- 
sipation region is 

(13) 



which corresponds to a thermal energy density in the 



3.3. Observational Constraints on the Jet Model 

Many of the basic physical features of the somewhat 
"canonical" jet model outlined above are strongly sup- 
ported and motivated by decades of observations. A syn- 
opsis of the corroborating evidence follows. 

3.3.1. Jet Kinetic Luminosity: L, and A^^^ 

Edge-on Fanaroff-Riley jets are better suited than 
head-on blazars to measure the jet kinetic power. For 
blazars, the highly Doppler-beamed, variable, broad- 
band radiation output must be properly modeled in or- 
der to extract physical parameters at the jet dissipation 
scale. At any given frequency, the value of the comoving 
radiation flux F^, itself is difficult to determine. In fact, 
in the observer frame. 
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where 6 is the Doppler-beaming factor (e.g., Lind fc| 
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jet axis. Clearly, due to the strong dependence on the 
Doppler-beaming factor S, an accurate measurement of 
the comoving radiation flux F^, , and therefore of intrin- 
sic jet properties such as and Aj^.^^ = LJL-^^^^, is 
difficult to perform. 

In the case of Fanaroff-Riley sources, the terminal 
sites of jet dissipation, the so called radio lobes, can be 
thought of as the thermal reservoirs of the initially or- 
dered and collimated jet mechanical energy. Since the 
energy injection rate outstrips the radiative cooling rate, 
radio lobes have no choice other than to expand. Plasma 
temperature and density within a radio lobe are typ- 
ically extrapolated from its X-ray continuum and line 
emission. From the inferred lobe enthalpy E^^^^, the jet 
kinetic power can be estimated as 



where dEf./dt is the rate at which a cold electron loses 
its orderly kinetic energy via the bulk-Comptonization 
process, cr^ is the Thomson cross section and JJblr = 
C -^acc /(^TT-R^c) is the energy density at radius R result- 
ing from the fraction ^ of accretion lumino sity L^^^ re- 
processed and isotropized by the BLR. Eq. (20 1 leads to 
an expression for the Thomson optical depth t^^j^ from 
the jet base up to the dissipation region: 
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rr'ri' , (21) 



E,._ 



t,. 



(18) 



where t^^^^^ ~ ^lobo/^s is the characteristic time required 
to inflate a lobe with width l,^^^^ at the sou nd speed Cs 



(e.g., Birzan et al.||2004 Shurkin et al. 
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At the powerful end, the lobes of J^'K 
average kinetic inputs of Lj ~ lO''^'^ 

sponding to the Thomson Eddington limit for black holes 
with mass 10^ — 10^° Mq.^ The fact that the upper limit 
for Lj is close to the Eddington limit for the largest black 
holes informs us that, at least during periods of power- 
ful activity, a value of A^.^^ ~ 1 may be appropriate for 



these sources, and that A„ 



1 may serve as a reason- 



able upper bound. For powerful FSRQs, thoug ht to be 
the on-axis c o unter part of massive FR II galaxies, Celotti| 
& Ghisellini (20081 infer comparable values for Lj from 
modeling the blazar spectral energy distribution. 

3.3.2. Optical Depth and Pair Content: t^.^^ and n^/up 

Information on particle composition at the jet dissipa- 
tion scale is best extracted from the spectra of FSRQs, 
where both jet and disk emissi on are present. T he gen- 
eral procedure is outlined by Sikora et al. (1997). 

The jet is bathed by optical/ U V photons from the BLR 
and IR photons from the presumed obscuring torus. Ow- 
ing to the large bulk Lorentz factor F ^ 10 of the jet flow, 
"cold" electrons in the jet have the capacity to Compton 
up-scatter these relatively soft photons (a process known 
as "bulk-Comptonization" or first-order kinetic Sunyaev 
- Zel'dovich effect) up to characteristic energies 



Vl keV 



(19) 



where hv,. 



10 eV is the typical seed photon energy. 



The expected bulk-Compton luminosity L-^^ from cold 
jet electrons may be written as a volume integral extend- 
ing from the jet base ~ Rc/^r^d the dissipation scale 
R^-^^, i.e., below the region where electrons are shock- 
heated: 



dVr. 



dE, 



dt 



2T^ / dVn.ca^T'^U^, 



(20) 



^ Black-hole mass is usually estimated from optical and UV lines, 
since their width supposedly measures the depth of the gravita- 
tional potential at the BLR. 



where we have imposed that the expected bulk-Compton 
luminosity L^^ should not exceed the observed inte- 
grated power Lgj, at soft X-ray energies (wh ere the bulk- 
Compton emission should peak, see eq. (19l). 



tgalaxies possess constramed to be 
10^°ergs~^, corre- Ue ^ 



Assuming that the flux of cold electrons is conserved 
along the jet (i.e., pair injection occurs primarily at the 
jet base) , it follows that the Thomson depth r^^^^ down to 
the dissipation scale satisfles t^.^^ < r^^j^, since the opti- 
cal depth is greater at smaller radii. Thus, from eqs. ( |12[ ) 
and (211 the electron to proton ratio of the jet flow is 
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(22) 



which informs us that electrons and positrons, though 
dominant by number, only advect with them an insignif- 
icant fraction of the jet power. That is, the kinetic power 
of t he l et is almost entir ely carried by ions, as assumed 
Also, from eq. (12) with rie/up 10 it follows 
0.2, and the now in the dissipation region is 



in ^3^ 
thaiT 
optically thin. 



3.3.3. Magnetic Field Strength: 

As in the case of estimating the jet composition, 
FSRQs are the best sources for measuring the relative 
strength of the magnetic fleld at the dissipation scale, 
parametrized by £3 — U'^/U[^. In the event that all of 
the low-energy synchrotron (with integrated luminosity 
Lg) and high-energy inverse Compton {L-^^) emission is 
powered by electrons accelerated in the dissipation re- 
gion, then 



U'. 



(23) 



where U'_^^^^ is the comoving energy density of seed pho- 
tons for the inverse Compton process. For FSRQs, U' . 
is mostly contributed by external photons from the BLfi 
and obscuring torus (as opposed to synchrotron seed pho- 
tons, that dominate in BL Lacs ) , so t hat U' ~ U, 



^^^^^^ , .^^ soft •-'BLR' 

where [/j3LR defined in J 3.3.2 With help from 



eq. ( 14 1 for the form of U^^ we have 
UL 2e, 



where we made liberal use of the fact that L 
our jet model in FSRQs. It follows that 



0.05 



0.1 



10- 



(24) 



for 



(25) 



which is consistent with the value e„ 



^ 0.1 chosen in 

^3.2 Comparable values result from detailed modeling 
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of blazar s pectral energy distributions (Celotti & Ghis- 
ellini|20d8 1 . This confirms that the electromagnetic con- 
tribution to the jet energy flux at the dissipation scale is 
negligible compared to the proton kinetic flux. 

The magnetic energy fraction computed in eq. ( 25 1 
is independent from the black hole mass M,. The ap- 
parent lack of dependence of on may be mis- 
leading. In fact, the fraction ^ of accretion luminosity 
that is re-processed and isotropized in the BLR will in 
principle depend upon M,, since the characteristic disk 

temperature scales as (x M »"^^^ at fixed Eddington ratio 
( Shakura fc Sunyaev|1973 1, which leads to a correspond- 
ing change in the photo-ionizing luminosity per unit of 
accretion power. Nevertheless, for the level of accuracy 
of this work we ignore such complications. 

3.4. Summary of the Jet Model 

Our "canonical" jet model may be considered more as a 
"consensus" jet model. That is, it is based upon decades 
of observations and modeling of radio-loud quasars, 
rather than any deep theoretical principle. We do not 
address the nature of the ultimate source of mechanical 
energy at the base of the jet. We simply exploit the fact 
that the jet is launched with some mixture of protons 
and electrons at a constant Lorentz factor F ^ 10 and 
half-opening angle 9 ^ 0.1, as often inferred in bright 
core-dominated sources. Although pairs may outnumber 
protons, the jet power is dominated by the ions. The 
most important feature of the jet model in terms of cos- 
mic ray production is that dissipation by internal shocks 
at a distance R^^^^ — r^^c/^iad from the black hole leads 
- almost certainly on theoretical grounds - to the ran- 
domization of bulk baryonic energy and - without doubt 
on observational grounds - of bulk leptonic energy. The 
shock-heated protons that escape the jet into the inter- 
stellar medium of the host galaxy may have profound 
implications for the galaxy evolution. 

The most important aspect of our jet model is that it 
is extremely simple with relatively few inputs. That is, 
we adopt typical values for F, AF, 9 and e^^^ and we de- 
rive physically-motivated estimates for A^,^^ , e^j^ and . 
From these parameters, we extract all the necessary in- 
formation required to model our jet-powered cosmic ray 
feedback mechanism, as we describe in the next section. 

4. DETERMINATION OF e^^ AND L^r FOR RADIO-LOUD 
QUASAR ACTIVITY 

In our model of quasar self-regulation, we rely upon 
interstellar cosmic rays generated in the radio core of 
quasar jets. As discussed in |3 2| in the jet dissipa- 
tion region a fraction e^j^ ^ 0.1 of the jet's bulk kinetic 
power is converted into random thermal form by internal 
shocks. The average comoving energy of shocked protons 
will be in the GeV range, since their mean Lorentz fac- 
tor is 7p ~ 1 -f e^,^ ; this corresponds in the host-galaxy 
frame to a kinetic energy ~ F 7^771^0^ ~ 10 GeV. We 
now address the question of whether a sufficient number 
of these moderately-relativistic cosmic rays can diffuse 
out of the jet into the host galaxy, thus providing the 
coupling between jet power and interstellar gas required 
for our feedback scenario. 

We find that most of the randomized protons produced 
by internal shocks are convected with the jet flow away 
from the dissipation region without suffering significant 



losses. Only a small fraction (~ t^dv/^dnt t) escape 
from the jet before being advected into regions where 
magnetic irregularities in the flow are too weak to allow 
for significant particle diffusion. Here, t_^^^ is the advec- 
tion time through the dissipation region and j is the 
diffusion time across the jet, as measured in the jet co- 
moving frame. ^ It follows that the fraction of jet kinetic 
luminosity L , ending up as interstellar cosmic ray power 
will be 



t 



(26) 



It is the purpose of this section to estimate the cosmic ray 
efficiency e^^ and the corresponding interstellar cosmic 
ray power L^^ = e^-^L-^. Comparison of L^-^^ with the 
Eddington limit in cosmic rays defined in eq. ^ will 
assess if enough cosmic rays are injected into the host 
galaxy to initiate a momentum-driven wind capable of 
removing the galactic gas, thus potentially self-regulating 
the black hole - galaxy co-evolution. 

4.1. Important Timescales for 1 — 10 GeV Cosmic Rays 
in the Jet Dissipation Region 

Table [Tlsummarizes the important comoving timescales 
for 1-10 CSIeV cosmic rays in the jet dissipation region; 
in boldface, the timescales most relevant for our model. 
As discussed below, we find that: i) randomized rela- 
tivistic protons are produced by internal shocks on a 
much shorter timescale than advection, which is itself 
the fastest loss process for 1-10 GeV ions; ii) cosmic 
ray diffusion out of the jet is a slow process compared 
to advection, implying that only a small fraction of the 
shock-accelerated protons will contribute to the interstel- 
lar cosmic ray population. 

4.1.1. Cosmtc Ray Production and Losses 

Although eq. implies that the dissipation region 
is coUisionless, the presence of even modest magnetic 
fields reasonably ensures that the transition from the fast 
cold pre-shock flow to the relatively-slow hot post-shock 
medium occurs within a few Larmo r scales (as observed 
for the Earth's bow shock by, e.g., |Tanaka et al.|p!983l 
Sckopke et al. 19901. In the jet comoving frame, this 
implies that internal shocks heat and randomize the in- 
coming cold protons on the short timescale 



t^ 



3.9 X 10 



-5 



F?, 



> ,1/2 1/2 1/2 



Af 



,.1/2 



(27) 

where = eB'/i^jprnpc) is the Larmor frequency of 
shock-heated protons with characteristic Lorentz factor 
7p; B' is the comoving magnetic field strength in eq. ( 16 1. 
We remark that our cosmic ray-driven feedback model 
does not require highly- relativistic protons, possibly ac- 
celerated at internal shocks by some Fermi mechanism. 
Rather, the bulk of the self-regulation effect results from 
"thermal" protons with comoving energies in the GeV 
range, produced in the dissipation region on the timescale 
^hcat computed above. For the sake of completeness, we 
also include in Table [T] the characteristic timescale t^^j^ 

^ We point out that all timescales in this section are computed 
in the jet comoving frame. 
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of diffusive shock acceleration.^ 

Production of ~ GeV cosmic rays on the timescale t^^^^ 
estimated above is much faster than advection through 
the dissipation region, which happens on a characteristic 
time 



t 



;4.9 X lO^e'l ^Tl9.iM,^g s ,(28) 



where we have made use of the fact that in the co- 
moving frame the dissipation region is roughly spheri- 
cal {R^.^^/T ^ R^._,^6). This is of the same order as the 
adiabatic cooling time due to the jet lateral expansion: 



t 



1 R. 



/3. 



4.9 X 10^ e-j ^ Ti s , (29) 



where f3^^^ ~ T9 is the jet expansion velocity in the co- 
moving frame. 

As Table [l] shows, advection and adiabatic cooling 
are the fastest loss processes for 1-10 GeV ions. The 
table also includes the characteristic timescales for the 
following processes: radiative cooling via synchrotron 
and inverse Compton emission, inelastic collisions with 
other protons {pp collisions, with cross section <Tpp ~ 
3 X 10~^^ cm'^ at energies of a few GeV and inelasticity 
kpp — 1/2) and inelastic collisions with background pho- 
tons (p7 collisions) - resulting in photo-meson produc- 
tion and/or pair production via the Bethe-Heitler effect 
(for an estimate o f the corresponding cross section, see 



Sikora et al. 



19871). We also discuss the limiting case in 
which the proton thermal energy stored in the dissipa- 
tion region (see U^^ in eq. ( 14 1 ) is efficiently transferred to 



the emitting electrons, wfllch then cool via synchrotron 
and inverse Compton at the observed bolometric lumi- 
nosity Lg^j ^ 10**^ ergs" the corresponding timescale 
is a conservative lower limit which holds regardless of 
uncertainties in the efficacy of energy exchange between 
protons and electrons. 

TABLE 1 

Important Timescales for 1 — 10 GeV Protons in the Jet 
Dissipation Region 



PHYSICAL PROCESS 


TIMESCALE (s) 


SYMBOL 


Shock Heating 


3.9 X 10"^ 7p 


*hoat 


Fermi Acceleration 


2.1 X 1027^=* 
6.9 X 10* 


(Kolm) 

'dsa 
■■dsa 


Advection 


4.9 X 10^ 


*adv 


Adiabatic Cooling 


4.9 X 10^ 


*ad 


Synchrotron and IC Cooling 


8.8 X 10^57-^ 




P — P Inelastic Collisions 


1.1 X 10" 




p — 7 Inelastic Collisions 


2.4 X 10^2 7-1 




p — > e Energy Transfer 


1.7 X 10** 




Diffusion 


1.1 X 109 7p"'^^ 
3.5 X 10** 


(Kolm) 

diff ..1 
(3tr) 

*dift,.l 



When evaluating the above timescales, we use A^.^^ 
= 0.1, € . = 0.1, e,. = 0.1, £„ = 0.1 and M. = 10 



: 1, r : 



10, 
We 



Diffusive shock acceleration occurs on a timescale t^. 



*adv/''^i' where Tj > 1 is the cosmic ray optical depth in the jet 
dissipation region (see j ]4.1.2| . Therefore, although not required in 
our model, there woulti be enough time to accelerate a fraction of 
the shock-heated protons to suprathermal energies before the flow 
is advected away from the shock. 



also assume that L^^/L^ 
relevant for our model. 



4.1.2. 



10. In boldface, the timescales most 



Cosmic Ray Diffusion out of the Jet 

In the tangled and inhomogeneous fields of the dissipa- 
tion region, a fraction of the shock-heated protons may 
efficiently scatter with resonant magnetic fluctuations 
and diffuse out of the jet before being advected away 
with the jet flow. The cosmic ray diffusion timescale 
across the jet is given by 



R.: 



*adv 



(30) 



where Tj is the cosmic ray optical depth. ^ The resonant 
magnetic fluctuations providing the cosmic ray scatter- 
ing may be embedded in the jet flow with, e.g., a Kol- 
mogorov spectrum, or generated in the dissipation re- 
gion by the accelerated cosmic rays themselves via the 
so ca l led "streaming instabifity" (e.g., Kulsrud & Pearce 
T969| |WentzeIlpJ74| . 

TTconcise review of both scattering processes and their 
ability to describe the cosmic r ay distribution function in 
the Milky Way is given in §3 of [Socrates et al.| ( |2008[ ). We 
assume that the scattering mechanisms which account for 
the interstellar cosmic ray optical depth in the Galaxy 
operate in the core of quasar jets as w e ll and we employ 
the same scalings as in Socrates et al. ( 2008 1 . This may 
seem as a huge extrapolation, but in tne absence of di- 
rect observational constraints on the cores of radio-loud 
quasars, we utilize this assumption for lack of a better 
choice. We now estimate the value of r, expected for the 
two scattering processes mentioned above. 

The cosmic ray optical depth in the presence of res- 
onant Kolmogorov turbulence at the Larmor scale of 
shock-heated protons (r^^ j = c/w^) may be written 



SBo 
B' 



R..6 



1/3 



Ao 



2/3 



(31) 



where Ao, (5i?o and B' are respectively the stirring scale 
of the magnetic turbulence, the fluctuation amplitude at 
the stirring scale and the jet magnetic fleld compute d in 
eq. (16 1. Assuming SB^) ^ B' and Aq ^ R^iss^^ *5q. (31 1 



gives 

(Kolm) 

T. : 



2.3xlO^;V3AV6^i/6^i^r-^/^M:/*^ , (32) 



and the corresponding diffusion time across the jet is 
t<""""c.l.l X 10%-i/3Ai/6ei/6 ,1/6 ,-1 Yl'^'e , mI'^ s . 

difi,.I 'P Edd B.-l th.-l rad,-l 1 "1 

In the case that the cosmic ray opacity is provided 
by magnetic fluctuations self-generated via the streaming 
instability, the optical depth may be written 

r;"' — c.7.1e-\/Vif , (34) 

J .11 B,-l th,-l ' ^ ' 



where , = B' / ^47rp^.^__ is the Alfven velocity in the 
dissipation region. The diffusion time across the jet is 

^ The hierarchy between the characteristic timescales for dif- 
fusion Jiff J ) , advection {t^^^ ) and diffusive shock acceleration 



(*dsa) is such that t^jff j - 
have t Jiff J >t^^^ >tusA- 



and since r, > 1 
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then 

(str) 
^ditt.J 



where we have made use of L, 



A 



3.5 X 10^ e-\/2 ^-i r? 6* , M. g s . (35) 

B,-l th,-l rad,-l -1 ^ ' 



By modehn g the spectral energy distribution of pow- 
erful blazars, Celotti & Ghisellini] p008 ) infer that the 
~ 0.5 — 10 GeV electrons responsible for the observed 
synchrotron and inverse Compton emission should be in- 
jected at internal shocks with a power-law energy spec- 
trum ?T.(7e) oc 7-" with a ~ 2.5. If the mechanism 
responsible for the electron acceleration to suprather- 
mal energies depends only on particle rigidity, as is 
the case for a Fermi process, then a comparable spec- 
tral index should describe the accelerated protons. For 
first-order Fermi acceleration at mildly-relativistic in- 
ternal shocks this implies, for ^ GeV protons, that 
T^/S.E/E ^ TjAF/F ~ 1, since the inferred particle spec- 
trum has nearly equal energy per logarithmic interval. 
Here, IS.E/E is the average fractional energy gain per 
acceleration cycle. It follows that ~ F/AF > 1. 

4.2. Cosmic Ray Efficiency and Luminosity L^^ 

Armed with a better knowledge of the physical pro- 
cesses relevant for 1-10 GeV protons produced in the jet 
core, we can now compute the fraction e^j, of jet power 
injected into the host galaxy as cosmic rays, and the cor- 
responding interstellar cosmic ray luminosity L^^ . 

Since the fraction of shock-heated protons that can 
diffuse out of the jet before advection takes place is 
~ *adv /''ditf J ^ l/'^'j; ^ shown in eq. (30 1, the cosmic ray 
efficiency in eq. ( 26 ) may be rewritten as e^j^ ~ e^j^ /Tj . 



If the cosmic ray scattering is due to resonant ma gnet ic 
fluctuations with a Kolmogorov spectrum (cf. eq. (32 1), 
then 



■-CR — ^(Kolm) 



4.4 X 



g5/6 pl/3 
IQ-S 1/3 th,-l 1 



,1/6 1/6 ^'^'.S 



1/6 



(36) 



where 7^ ~ 1 — 10 may be thought of as the characteris- 
tic comoving Lorentz factor of shock-accelerated protons. 
Instead, if the cosmic ray opacity is provided by mag- 
netic waves self-generated via the streaming instability 
(cf. eq. (34)), the resulting cosmic ray efficiency is 



1.4 X 10-2 £1/2 g3/2 

B,-l th,-l 



(37) 



As eqs. (36 1 and (37 1 suggest, the cosmic ray coupling 
efficiency e^j, = L^^L^ is either weakly dependent (for 
Kolmogorov turbulence) or not dependent at all (for the 
streaming instability) on the black hole mass M,. In 
other words, our cosmic ray-driven feedback scenario is 
manifestly self-similar or scale-independent per unit of 
jet power (or equivalently, per unit of black hole mass, 
for fixed radiative efficiency). 

The cosmic ray luminosity L^^^ — Ccr^i injected into 
the host galaxy is then, for the two scattering mecha- 
nisms discussed above. 



C""-5.7xlO-7y 



A5/6 5/6 pl/3 

1/3 Edd th,-l 1 



M.%%rgs-i ,(38) 



1/6 



C^l.8 X 10^^ A_eV2 e3/2^M.,9 ergs-^ , (39) 



Edd Edd, • " 



Comparison of eq. ( 39 1 with the Eddington limit 
in cosmic rays defined m eq. ^ (or equivalently, of 
eq. (37 1 with eq. (l7|) suggests that, if the resonant 
magnetic turbulence for cosmic ray scattering in the jet 
core is mainly generated by the streaming instability, 
there may be enough power in interstellar cosmic rays 
to launch a momentum-driven wind which could eject 
the galactic gas. Although the momentum requirement 
^CR ^ ^Edd CR is ^ necessary condition for efficient self- 
regulation of the black hole ~ galaxy system, it is not 
sufficient by itself; in fact, we must also require that the 
time-integrated cosmic ray energy injected into the host 
galaxy during this epoch of super-Eddington cosmic ray 
activity should be comparable to the binding energy of 
the gas. In f|5]we address this important issue. 

4.3. Cosmic Ray Propagation within the Host Calaxy 

In the frame of the host galaxy, the shock-accelerated 
protons which diffuse out of the jet are relativistically 
beamed along the jet axis within an angle ~ l/F. Since 
the cosmic ray optical depth for scattering off resonant 
magnetic fluctuations in the galaxy's interstellar medium 
is very large (r^ ~ 10^, as we show below), the cos- 
mic ray momentum distribution quickly isotropizes so 
that their overall effect on the galactic gas resembles a 
spherically-symmetric pressure force in the direction op- 
posite to the galaxy's gravitational center. As the in- 
terstellar cosmic rays diffuse towards larger scales, their 
energy does not appreciably change from the value ~ 
Tjpnipc'^ ~ 1 — 100 GeV of their birth, since the average 
fractional energy loss per scattering is only ^ t^"^ <^ 1 



(e.g., Kulsrud fc Pearce 1969 Wentzel 1974). Instead, 



the fractional momentum change per scattering is ~ r" , 
which allows cosmic rays to be extremely effective in pow- 
ering a momentum-driven outflow of interstellar gas. In 
this section, we study the cosmic ray propagation within 
the host galaxy, providing an estimate for the cosmic ray 
optical depth Tg and diffusion time t^.^j ^ . 

The characteristic galactic scale radius Rg for an 
isothermal sphere is 



Rn — 



at 



0.76 ( 



*,300 



'*,300 



kpc , 
(40) 

if the black hole - galaxy system follows the M, — Mi, 
and Af, — ct^, relations. We take Rg as the characteris- 
tic distance where cosmic rays produced in the jet core 
interact with the galactic gas. 

If we pin the ratio between cosmic ray mean free path 
and Thomson mean free path at the galactic scale Rg to 
the value A^^^/A^ ~ 10-^ appropriate for the Milky Way, 
the cosmic ray optical depth up to Rg is 



1.5 X 103/,,.i ( ) , (41) 



where pg is the gas mass density at Rg for an isothermal 
sphere. The cosmic ray diffusion time up to Rg is then 



^dm,g — '^a 



8.9 X 107.-1 
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We expect that, during epochs of vigorous jet-powered 
cosmic ray activity, the galactic structure will be appre- 
ciably affected by the cosmic ray pressure on timescales 
compa rabl e to the cosmic ray diffusion time computed 
in eq. (42 1, ft follows that, if the radio jet is powered 



by gas accretion onto the central black hole, then i^.^j ^ 
might be a reasonable upper limit for the duration of the 
radio-loud phase. In j]5.2| we further comment on this 



and examine the implications of the scaling t^.^^ ^ oc 



As discussed by Socrates et al. (20081, the cosmic ray 



pressure force may be reduced by ^/ upp^g 



/tditi.g if the dif- 
fusion timescale t^^^^ g is significantly longer than the time 
tpp^g required to deplete the cosmic ray energy via inelas- 
tic scattering with background protons resulting in pion 



pp 



production. For a cross section cr, 

GeV energies and a scattering inelasticity Kpp 
we find 



1/2, 



nip 1 

KppCTppC Pg 



1.8xl08/-.\<30oyr, (43) 



which is marginally longer than t^. j,. ^ , so that we can 
ignore losses due to pion production.^ 

Finally, we show that cosmic ray optical depths compa- 
rable to the value in eq. (41) may be derived by making 
the extreme assumption mat the magnetic energy den- 
sity in the galaxy is roughly in equipartition with the gas 
random kinetic energy density, i.e. i?^/87r ~ 3/2 PgCr^ at 
the characteristic radius Rg. If the cosmic ray scattering 
is provided by background magnetic turbulence with a 
Kolmogorov spectrum, the optical depth up to Rg is 



^^0 

B 



Rn 



1/3 



Ro 



2/3 



4.4xio3/;^/^S„ 



(44) 

where Aq ~ Rg is the stirring scale and 5Bq ~ B the 
fluctuation amphtude at the stirring scale for the inter- 
stellar magnetic turbulence, and r-^ g is the Larmor ra- 
dius in the galactic magnetic field B for a, ^ 10 GeV 
proton. Instead, if the cosmic ray opacity results from 
Alfven waves self-generated via the streaming instability, 
the cosmic ray optical depth will be 



5.8 X W'a-l,, 



(45) 



where ^ = B / yjAnpg is the Alfven velocity in the 
galaxy. Note that both scattering mechanis ms y ield an 
optical depth comparable to the result in eq. (41 1. 



5. COSMIC RAY FEEDBACK RESULTING FROM 
EPISODIC RADIO-LOUD ACTIVITY: ORIGIN OF THE 
M. - A-h RELATION 



As discussed in §4.2| during epochs of powerful jet ac- 
tivity the cosmic ray luminosity L^^^^ injected into the 
interstellar medium of the host galaxy may exceed the 
galaxy's Eddington limit in cosmic rays L-^^^ ^-^ (com- 
pare eqs. (39 1 and (H). Due to the outward cosmic 
ray pressurerorce, hyorostatic balance is lost and a cos- 
mic ray momentum-driven wind develops, removing the 
fuel for further star formation and black hole accretion. 
However, the bulge and black hole growth is completely 
choked off only if the radio- loud phase lasts long enough 
such that the total energy output in cosmic rays, which 
would eventually couple to the interstellar medium, is 
comparable to the binding energy of the galactic gaseous 
component. In other words, the momentum require- 
ment L^pf > Lj,jj is a necessary prerequisite, but ef- 
ficient self-regulation of the black hole - galaxy system 
is achieved only when the energy requirement discussed 
above is fulfilled as well. 

A generic feedback mechanism is capable of self- 
regulating the combined black hole - galaxy growth only 
if the total energy injected into the interstellar medium 
Ai?j^j is comparable to the gravitational energy of the 
gas: 



AE. 



Eg ~ fgNhat 



(46) 



where we have adopted an isothermal sphere. In what 
follows, we assume that the gas fraction fg is indepen- 
dent of black hole mass and stellar velocity dispersion. 
We now contrast our jet-powered cosmic ray feedback 
model with self-regulation mechanisms that act during 
the optically- luminous "quasar phase." 



5.1. Black Hole Self- Regulation during the Luminous 
Quasar Phase 



The Soltan ( 1982 ) argument, along with the work of 
Yu & iremaine '( 2002 1 , indicates that supermassive black 
holes at the center of galaxies build up their mass primar- 
ily by an act of radiatively-efficient accretion, in a rela- 
tively short-lived high-luminosity "quasar phase" . Dur- 
ing this epoch, the accretion luminosity approaches the 
black hole Thomson Eddington limit L^^^ 

If energy release during the optically-bright quasar 
phase is responsible for black hole self-regulation, the to- 
tal energy Ai?Qp absorbed by the galactic gas during the 
quasar lifetime Ai„p should satisfy AiJ„p ^ Eg, as pre- 



scribed by eq. (|46| . If only a fraction e^p of the accretion 



luminosity L_^^^ can couple to the interstellar medium of 



^ We point out that cosmic ray losses due to pp collisions are 
negligible all the way down to pc scales. If the Thomson op- 
tical depth in the BLR is ~ 0.1, the extreme assumption 
that the cosmic ray diffusion time in the BLR be comparable or 
longer than the proton-proton collision time implies a CR opti- 
cal depth at pc scales > mpKes/C^ppCTpp) ~ 450. In 
the BLR, the gravitational force is still primarily provided by the 
black hole, and the correspoding Eddington luminosity in CRs will 
be iEdd.CR = 'f^Edd,.^T/T-CR - 2.0 X 10«Af.,g ergs"!. This is 
smaller than the CR luminosity in eq. l |39| . Therefore, in this 
case, the CR flux would disrupt the entire tSL/R. Since the BLR ex- 
ists, we can comfortably reject the extreme assumption that CRs 
are significantly destroyed at pc scales. 



gu- 
ianl 



^ Caution must be taken when applying thelSoltan'('l982 ) ar; 
ment to black hole demographics, encapsulatecT by the Magorr: 
|et al.| |T998j l relation. The luminosity and mass function of super 
massive black holes peak close to a black hole mass M, ~ 10* Mq 
(see dotted line in Fig.jlJ. Relatively small black holes, like the one 
at the center of the Milky Way, are both too rare and faint, such 
that a Soltan-like argument cannot be employed in order to deter- 
mine whether or not their mass was built up by radiatively-efficient 
Eddington-limited accretion. However, it seems likely that rela- 
tively small black holes with masses as low as M, ~ lO'' - 10'' Mq 
can radiate close to the Eddington limit, which therefore suggests 
that the build-up of black hole mass in relatively small systems 
takes place during a bright short-lived "quasar phase" as well. 
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the host galaxy, we require 



^^QP - ^QP ^^QP 



Assuming L,^^^ ~ -^Edd . ; this can be rewritten as 



M, 



in G rup c 



^QP ^'•QP 



(47) 



(48) 



where Ai?j ~ Aij^pLj is the time-integrated kinetic out- 
put of the radio jet, At^p is the duration of the radio-loud 
phase and ej^p is the efficiency for coupling the jet power 



Lj to the galactic gas. The energy balance AiJ^^p 
in eq. ( 46 ) requires that 



The quasar hfetime At^g , set by the Soltan ( 1982 ) argu- 
ment, is roughly comparable to the Salpeter time of the 
black hole: 



^^RP - CrP ^^RP ^Edd^Edd,. ^ Eg ' 

where Aj^.^^ = L-^jL^^^ This can be rewritten as 



M, 



47r G nip c 



^RP ^''RP ■'^Edd 



En 



(52) 



(53) 



AE, 



47r G rrir 



i4.5 X lO^e^^^ ^ yr ,(49) from which it is apparent that the product e^p Aij^p A 



where AE, ~ e^_^^M,c^ is the total radiative energy out- 
put of the accretion process. Since e^^^ depends only 
on black hole spin, the quasar lifetime At^^ ~ ig^,p 
should be independent of black hole mass. It follows 
from eq. (48) that, in order for a system to lie on the 



Af, — Mi, relation, the coupling efficiency Cqp in the lu- 
minous quasar phase must depend on the stellar velocity 
dispersion ct* such that 

Eqp cx al . (50) 

This implies that, if the quasar phase were to self- 
regulate the black hole and galaxy evolution, the feed- 
back mechanism at work would not be a universal, scale- 
free or self-similar process. We stress again that this 
follows from Ai^p ~ ^saipi or equivalently from the fact 
that the total quasar energy output At^p L^^^ is observa- 
tionally pinned to be close to the maximum amount of 
energy released in the formation of the black hole, given 
by AE,.^ 

There is no apparent reason why the coupling efficiency 
Eqp may not depend on the stellar velocity dispersion 
(T^ or the black hole mass M,. However, it should de- 
pend upon a^, and M, in such a way that makes the 
ratio M,/AIi, a constant across nearly four decades in 
black hole mass. This is, of course, not impossible, but 
it would have to involve a cosmic conspiracy with re- 
spect to the gas dynamics of black hole self-regulation. 
We now show how our jet-powered cosmic ray feedback 
scenario provides a satisfactory solution for this apparent 
contradiction. 

5.2. Black Hole Self- Regulation During the Explosive 
Radio-Loud Phase 

For a generic feedback mechanism operating during 
epochs of AGN radio-loud activity, the total energy 
AE-p^p injected into the interstellar medium is 



AE^ 



.AE, 



Ej^p Aip^pZ/j 



(51) 



° The scaling in eq. seems to be in contradiction with 

the simulations by Di Mattco ct al. (2005), which reproduce the 
M, — Oi, relation assuming that a constant fraction of black hole 
accre tion energy is deposited in the galacti c ga s i jSpringel et al.| 
20051. However, we remark that e^p in eq. l|50| is tne traction ot 



accretion energy available to unbind the galactic gas, and not just 
the fraction of energy deposited in the gas. If the simulations by 
Di Matteo et al. (2005) are taken at face value, then there is a hid- 
den parameter that limits the fraction of deposited energy which is 
available to unbind the galactic gas. The hidden parameter should 
scale as oc ct? . 



for radio-loud phases replaces e^p Ai^p in quasar epochs 
(cf. eq. (48 1). For a black hole - galaxy system that lies 
on the AiT"— M^, relation, this implies 



^RP ^^R.P -^Edd 



(54) 



As o pposite to the quasar feedback scenario discussed in 
fj5.1| eq. (54 1 does not immediately constrain the cou- 



pling efficiency ej^p to depend upon the stellar velocity 
dispersion. The possibility that the feedback efficiency 
gpp is independent of M, and cr* would then require that 
Atpp ^EdA increases with the stellar velocity dispersion, 
or that the time-integrated jet kinetic output per unit of 
black hole mass should scale as oc trj. 

If black hole self-regulation in the radio-loud phase is 
mediated by cosmic rays produced in the jet core, we 
have shown that the coupling efficiency e^^j^ is indeed in- 
dependent (for cosmic rays scattering with Alfven waves 
self-generated via the streaming instability, see eq. (37)) 
or weakly dependent (for resonant diffusion by magn etic 
turbulence with a Kolmogorov spectrum, see eq. (36 1 ) on 
the black hole mass.^ In our cosmic ray-driven feedback 
model we should then expect that At^^^ -^Edd ^ '^^ or, if 
Ajjj^j is a constant with black hole mass, that the dura- 
tion of the radio-l oud phase should scale as Atj^p c>c . 

As discussed in j ]4.3| a reasonable upper limit for At^^ 
may be given by the cosmic ray diffusion timescale t^^^^ ^ 
within the host galaxy. In fact, if the jet is powered by 
black hole accretion, AGN radio activity would be termi- 
nated when the interstellar cosmic rays have diffused up 
to the scales where most of the gas resides, and started to 
push it outward, thus preventing further accretion. In- 
terestingly, if fg and A(,j^/A^ do not depend on black hole 



mass, eq. (42) shows that t^.^^g cx a^, the same scaling 
that At^p sEould have in order to satisfy eq. ( 54 ) with a 
scale-independent coupling efficiency. 
Note that for At^^ ~ t ^.^^ ^ , eq. ( 42 1 suggests that 



the duration of the radio-loud phase for the most mas- 
sive galaxies may be comparable to the Salpeter time in 
eq. ( 49 1 , whereas in relatively small galaxies it should last 
less than the quasar phase. In other words, for low-mass 
galaxies the black hole self-regulation would be confined 
to an episodic radio-loud epoch much shorter than the 
time required to accrue the black hole mass. 

" We have defined as the fraction of jet power injected into 
the host galaxy as cosmic rays. However, due to th e large cosmic 
ray optical depth within the host galaxy (see i ]4.3| , the whole in- 
terstellar cosmic ray power will eventually be transferred to the 
galactic gas, so that t^,p is also a good proxy for the feedback 
coupling efficiency. 
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Fig. 1. — Black-hole mass function (M, (f>{M,); dotted line) and 
kinetic energy output of radio jets integrated over cosmic time 
(Ai?jAf, (^{M,); solid line), versus black hole mass Af,. The func- 
tion (f){M,), namely the number density of black holes with mass 
in the interval [M, , M, -l-dAf,], is derived from the stellar velocity 
dispersion function of early-type galaxies by Shcth ct al. (2003) 
via the M, — c* relation. The jet time-integrated kmetic output 
AS J is supposed to satisfy the energy balance tj^p Ai?j ~ fgMi^a^ 
requir ed fo r black hole self-regulation during radio-loud phases, as 
in cq. I52I1. We have assumed a coupling efficiency e^p = 10"'^, as 
found tor our cosmic ray-driv en fe edback rnodel (somewhat an in- 
termediate value between eq. |36| and eq. l |37| ) and a gas fraction 
fg = 0.1 for the host galaxy; the black hole - galaxy system lies on 
the M, — Af* and Mm — ct* relations. We set hyo to be the Hubble 
constant Hq in units of 70 kms~^ Mpc~^. 

6. DISCUSSION AND COMPARISON WITH OBSERVATIONS 

We propose a self-regulation mechanism for supermas- 
sive black holes and their host bulges that operates in 
actively-accreting systems with powerful radio-loud ac- 
tivity. In the core of relativistic radio jets, internal shocks 
resulting from the dissipative interaction of shells of mat- 
ter, intermittently ejected from the central engine, con- 
vert a fraction of the ordered kinetic energy of the jet 
flow into thermal form. If a sufficient number of the 
protons (or "cosmic rays" ) heated and randomized at in- 
ternal shocks escape the radio core into the interstellar 
medium of the host galaxy, they may profoundly affect 
the evolution of the galaxy and its black hole. To quan- 
tify their effect on the hydrostatic balance of the galac- 
tic gas we define an Eddington limit in cosmic rays for 
the host galaxy. For a phenomenologically-motivated jet 
model, we show that during powerful radio-loud phases 
the power in 1-100 GeV interstellar cosmic rays is large 
enough to break the cosmic ray Eddington limit for the 
host galaxy, so that momentum balance of the galac- 
tic gaseous component is lost and a cosmic ray-driven 
wind develops, that removes the galactic gas. If this 
super-Eddington cosmic ray activity lasts long enough, 
the time-integrated cosmic ray energy input into the in- 
terstellar medium may exceed the binding energy of the 
gas, and the whole galaxy's gaseous phase will become 



unbound. In doing so, any fuel for further black hole 
growth and star formation will be removed, thus affecting 
the combined evolution of the black hole and its stellar 
b ulge, as implied b y the M « — M^, and M, — cr^, relations. 



Morganti et al. ( 2005 1 report the detection of fast 



(~ 1000 kms~^) large-scale (~ 1 — 10 kpc) massive out- 
flows of neutral hydrogen in several powerful radio galax- 
ies. They claim that the outflows are spatially associated 
with radio knots extended along the jet, but this may just 
result from the fact that H I absorption can be traced 
only in the presence of a bright background radio contin- 
uum. Fast large-scale outflows of low-ionization species 
with nearly spherical morphology have been d etected in 
the powerfu l radio source MRC 1138-262 by Nesvadba] 



et al. (|2006| . The fact that the outflow is "cold"^ (atomic 
or weakly ionized) and "fast" (super-virial) implies that 
it is the result of catastrophic momentum, rather than 
energy, exchange with the interstellar medium. In other 
words, it is possible that some sort of Eddington limit 
is being brok en. Among the moment um-driven feedback 



mechanisms, Nesvadba et al. ( 2006[ ) concl ude that nei 



ther radiat i on-powered quasar w inds (e.g., Fabian|[T999 



King 2003 Murray et al. 20051 nor direct coupling be 



tween the jet and the interstellar medium in the "den 



tist drill" model (Begelman & Ciofli 19891 seem suffi- 
cient to explain the large-scale gas kinematics observed 
in MRC 1138-262. Direct coupling between the jet and 
the galactic gas should result in the outflow being con- 
fined within a narrow cone around the jet axis, contrary 
to observations. Regarding the effect of radiation pres- 
sure from t he quasar photon ou tput on interstellar dust 
grains (e.g., Murray et aH2005 1, it is hard to understand 
how this could remove a significant fraction of the galac- 
tic gas, since the UV photons will be degraded to IR 
frequencies after a few scatterings, within a parsec from 
the black hole. Also, the optical depth for the resulting 
far IR photons will not exceed unity beyond a few tens 
of parsecs. 

The jet-powered cosmic ray feedback scenario pre- 
sented here does not suffer any of these deficiencies. The 
cosmic ray luminosity injected into the host galaxy is 
smaller than the black hole photon power by roughly 
two orders of magnitude, but cosmic rays exchange mo- 
mentum with the galactic gas ~ Tg/r^^ ~ 10^ more ef- 
ficiently. Here, r^^ ~ 1 is the optical depth for UV 
photons on dust grains. Moreover, since the cosmic ray 
optical depth in the galaxy's interstellar medium is very 
large {tq ~ 10"^), any memory of the cosmic ray momen- 
tum distribution at injection is quickly lost and their 
effect on the galactic gas should resemble a spherically- 
symmetric outward pressure force, so that the resulting 
gas outflow is nearly spherical. Also, the fraction of cos- 
mic ray energy lost at each interaction with the inter- 
stellar gas is minimal (~ t"^), meaning that they can 
propagate up to the large scales where most of the gas 
resides without suffering signiflcant losses. 

An interesting result of our study is that, per unit of 
black hole energy release, the cosmic ray feedback effi- 
ciency e^^j^, which gives the fraction of jet power available 
to unbind the galactic gas, is roughly a constant with 
black hole mass. Instead, for self-regulation models rely- 
ing on the black hole photon output d uring radiatively- 
efficient quasar phases (e.g.. Silk fc Rees||1998| [Fabian] 
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T999|ICiotti fc Ostriker|200Tl|King|2003|IDi Matteo et al 



2005||Murray et al.||2005| [Hopkins et al.t2006p , the feed- 
back efficiency is constrained to scale as ^ Eg/ AE, oc ctJ 
in order for the system to he on the M, — M^, rela- 
tion. Here, Eg ~ fgMi^al is the binding energy of the 
galactic gas and Ai?, is the time-integrated energy out- 
put resulting from black hole accretion. It would be a 
surprising coincidence if an intrinsically scale- dependent 
self-regulation mechanism were to result in the scale-free 
M, — relation, which holds for nearly four decades in 
mass. 

The explosive radio-loud phase does not suffer from 
such severe constraints. Rather, since the coupling ef- 
ficiency e^j, for our cosmic ray-driven feedback scenario 
is scale-independent, the energy balance e^^j^AiJj ~ Eg 
required for black hole self-regulation implies that the 
time-integrated jet kinetic output should scale as Ai?j oc 
Mi,a1 for a constant gas fraction fg. By inferring the 
black hole mass function M, (j){M,) (dotted line in Fig.[l]) 
from t he velocity di spersion function of early-type galax- 
ies by Sheth et al. (2003) via the M, — cr^, relation, we 
can predict the dependence on black hole mass of the ki- 
netic energy output of radio jets integrated over cosmic 
time (A£^jM, (/)(M,); solid line in Fig. [ij, for systems 
lying on the M, — M^, and M, — ct^ relations. We find 
that its slope at the low-mass end is ~ 2.5 and it peaks 
at M, ~ 3 X 10^ M0. With more reliable measurements 
of jet kinetic power and black hole mass, this could pro- 
vide a stringent observational test for our proposed self- 
regulation scenario. For example, finding a slope close 
to 2.5 at the low-mass end would imply that the ratio 
between the time-integrated jet kinetic output and the 



binding energy of the galactic gas should be independent 
of black hole mass, thus strongly implicating that black 
hole self-regulation occurs in the radio-loud phase, irre- 
spective of the actual coupling mechanism itself. 

For a generic feedback mechanism acting during radio- 
loud epochs, the scaling AiJj (x M^trJ can be recast as 
Aij^pAj.^^ (X al, for a black hole - galaxy system which 
follows the M, — Mi, relation. Here, Aj^.^^ is the ratio 
of the jet kinetic power to the black hole Thomson Ed- 
dington limit and Aij,p is the duration of the radio-loud 
phase. This implies that, for systems where a significant 
amount of mass is being built up, radio-loud signatures 
will most likely be observed in galaxies with large black 
holes and bulges, in line with the b asic phenomenology of 
radio-loud AGN and quasars (e.g., Laor|200^ Best et al. 
2005 1 . The Af, — Mi, relation perhaps is then connected 
to the fact that radio-loud actively-accretin g objects are 
relati vely absent in relatively small systems ( [Greene et al. 
20061: black hole self-regulation takes place during the 
radio-loud phase and, due to the abundance of total en- 
ergy available, this phase is short-lived in relatively small 
systems because the amount of energy required to unbind 
the galaxy is relatively low. 
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